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PHASE FIELD MODELS FOR THIN ELASTIC STRUCTURES WITH 
TOPOLOGICAL CONSTRAINT 

PATRICK W. DONDL, ANTOINE LEMENANT, AND STEPHAN WOJTOWYTSCH 


Abstract. This article is concerned with the problem of minimising the Willmore energy in 
the class of connected surfaces with prescribed area which are confined to a small container. We 
propose a phase field approximation based on De Giorgi’s diffuse Willmore functional to this 
variational problem. Our main contribution is a penalisation term which ensures connectedness 
in the sharp interface limit. The penalisation of disconnectedness is based on a geodesic distance 
chosen to be small between two points that lie on the same connected component of the transition 
layer of the phase field. 

We prove that in two dimensions, sequences of phase fields with uniformly bounded diffuse 
Willmore energy and diffuse area converge uniformly to the zeros of a double-well potential 
away from the support of a limiting measure. In three dimensions, we show that they converge 
Ti 1 -almost everywhere on curves. This enables us to show T-convergence to a sharp interface 
problem that only allows for connected structures. The results also imply Hausdorff convergence 
of the level sets in two dimensions and a similar result in three dimensions. 

We furthermore present numerical evidence of the effectiveness of our model. The implemen¬ 
tation relies on a coupling of Dijkstra’s algorithm in order to compute the topological penalty 
to a finite element approach for the Willmore term. 
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1. Introduction 

In this article, we consider a diffuse interface approximation of a variational problem arising 
in the study of thin elastic structures. Our particular question is motivated by the problem of 
predicting the shape of certain biological objects, such as mitochondria, which consist of an elastic 
lipid bilayer and are confined by an additional outer mitochondrial membrane of significantly 
smaller surface area. 

In order to describe the the locally optimal shape of thin elastic structures, we will rely on 
suitable bending energies. A well-known example of such a variational characterisation of biomem¬ 
branes is given by the Helfrich functional IHel731 ICan70l 

( 1 - 1 ) £ h (£)= [ XH(H-H 0 ) 2 + x K KdH 2 

Jy. 

where £ denotes the two-dimensional membrane surface in R 3 and H and I\ denote its mean and 
Gaussian curvatures. The parameters xh, Xk and Hq are the bending moduli and the spontaneous 
curvature of the membrane and are, in the simplest case, constants. The integral is performed 
with respect to the two-dimensional Hausdorff measure Ti 2 . The special case when Hq,xk = 0 
and xh = 1 is known as Willmore’s energy 

(1.2) W(£) = [ H 2 dn 2 . 

J £ 

While our results extend to certain Helfrich-type functionals of more general form (see section [5]), 
the basic task we have set for ourselves is the following: Minimise Willmore’s energy in the class of 
connected C 2 -surfaces which are embedded into a domain such that the embedding has prescribed 
surface area S > 0. This problem does not necessarily admit a solution, but a minimiser exists in 
the varifold closure of these surfaces [ DMRIT . 

In this version of the problem, the outer mitochondrial membrane is the rigid boundary 912, but 
the results also apply if an elastic container is used, see section [5] 

We propose a phase field approach to this problem. For computations, this has some advantages 
to the naive approach of taking a gradient flow of Willmore’s energy (or an area preserving version 
of it). 

First, Willmore flow is a fully non-linear degenerate parabolic PDE of fourth order on a moving 
surface, which makes both analytic and numerical treatment difficult. In contrast, our phase- 
field flow equations are fourth order quasi-lincar parabolic with constant coefficients for the fourth 
derivatives. The price we have to pay for this relative luxury is, as usual, that we have to solve 
equations in one higher dimension. 

Second, Willmore flow is a fourth order PDE and as such does not admit a maximum principle. 
In particular, a surface originally confined to even a convex container could easily leave it along 
the gradient flow, whereas for us the confinement to Q is simply given by the choice of boundary 
condition and the domain of the phase field. 

The main challenge for phase fields is to control the topology of a limiting surface (in some sense) 
from quantities of diffuse interfaces. In this paper, we give a suitable concept of connectedness 
that admits a rigorous variational statement as £ -A 0 and can be implemented efficiently in a 
simulation. 

Heuristically, we introduce a quantitative notion of path-connectedness for the interfacial region 
of a phase field. The disconnectedness for two points x,y in the interface layer is measured by 
a geodesic distance (the infimum length of connecting paths) where walking in the interfacial 
region does not add any length, but leaving the interface adds a positive amount. The total 
disconnectedness is then computed as a double integral over the distances of interface points. 
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This functional poses analytic challenges since it is not a priori clear how ’two-dimensional’ the 
interfacial layer of a phase field is. Since the distance function uses the length of curves, away from 
a limiting surface we need the phase-fields to converge almost everywhere on curves, i.e. on objects 
of codimension 2. Also, the structure of the interface has to be understood in a precise way in 
order to estimate the aforementioned double integral. More rigorous statements are given below. 

Our analysis builds on ideas in IRS06 . Along the proof of our Theorem, we obtain a number of 
useful technical results about the convergence of phase field approximations to the limit problem 
which appear to be new in this context. Namely, for sequences along of uniformly bounded diffuse 
energy, in two dimensions we prove convergence of the transition layers to the support of a limiting 
measure in the Hausdorff distance and uniform convergence of the phase field to ±1 away from the 
limit curve. This is reminiscent of a result for minimisers of the Modica-Mortola functional among 
functions with prescribed integral [CC95I and more generally for stationary states 1 l'l 1)0 . 

In three dimensions, we show that neither result is true, but prove id-convergence on curves and 
that spt(/z) is contained in the Hausdorff limit of interfacial regions. This enables us to precisely 
control our topological energy term. We also show that phase fields are uniformly L°°-bounded in 
terms of their Willmore energy, also in three dimensions. 

The numerical implementation of our functional will be discussed further in a forth-coming 
article [DW161 . We use a variant of Dijkstra’s algorithm similar to the one of jBCPSlO] to compute 
the geodesic distance function used in the topological term of our energy functional. The weight 
in the geodesic distance is chosen to be exactly zero along connected components of the transition 
layer, which implies that it only needs to be computed once per detected connected component. 
This makes the functional efficient from an implementation point of view. 

The article is organised as follows. In the remainder of this section we give some further context 
for our problem. In section [2] we give a brief introduction to phase fields for Willmore’s problem 
in general (sections 2.1 and 2.2) and to our approach to connectedness for the transition layers 


(section 2.3). Our main results are subsequently listed in section 2.4 

Section [3] is entirely dedicated to the proofs of our main results. We directly proceed to show 
T-convergence of our functionals (section 3.1). In sections 3.2 and 3.3 we produce all the auxiliary 
estimates that will be needed in the later proofs. These are then used to show convergence of the 
phase fields, Hausdorff-convergence of the transition layers, and connectedness of the support of 
the limit measure (section [3.4[ ). 

In section [4] we show numerical evidence of the effectiveness and efficiency of our approach. To 
that end, we compare it to diffuse Willmore flow in two dimensions without a topological term 
and to the penalisation proposed in l )MH 1 1 . We have been unable to implement the functional 
developed in |DMR14] in practice, so a comparison with that could not be drawn. In section [5] we 
discuss a few easy extensions of our main results. 


1.1. Topology and Phase Fields. An often cited advantage of phase fields is that they are 
capable of changing their topology; in that sense our endeavour is non-standard. It should be 
noted that our phase fields may still change their topology (at least in three dimensions), only 
connectedness is enforced. 

Examples of topological changes and loss of connectedness in simulations for biological problems 
governed by bending energies or our type are given in [DW071 IDulOl . 

In BLS15] . a geodesic distance function has been used to minimise the length of a connected 
set K containing a prescribed set of points x\,... ,Xn in two dimensions (Steiner’s problem). Our 
setting is different in two ways: 1. Steiner’s problem has a finite number of a priori known points 
which need to be contained in K while the transition layer of the phase field has no special points 
and 2. the phase field approximation of Steiner’s problem works in dimension n = 2, while we 
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work in ambient space of dimension n = 2,3 where the curves used in the definition of the distance 
function have codimension 2. 

Previous work in |DMR11] provides a first attempt at an implementation of a topological con¬ 
straint in a phase field model for elastic strings modelled by the one-dimensional version of the 
Willmore energy, Euler’s elastica. This technique prevents transitions in simulations for simple 
situations, but may fail in more complex cases, see section [4j Similar numerical approaches to 
controlling the topology using a diffuse Euler number are discussed in (DLW051IDLRW07] . 

This approach was complemented by a method put forth in [DMR.14] . which relies on a second 
phase field subject to an auxiliary minimisation problem used to identify connected components 
of the transition layer. While our functional can be seen as using a diffuse measure of path- 
connectedness, the functional in [DMR.14j generalises more directly the notion of connectedness. 
For this model, a T-convergence result was obtained, showing that limits of bounded-energy se¬ 
quences must describe a connected structure. Unfortunately, the complicated nested minimisation 
problem makes it unsuitable for computation. 

Approaches of regularising limit interfaces have been developed by Bellettini in [Bel97l and 
investigated analytically and numerically in |ERR14j . The approaches work by introducing non¬ 
linear terms of the phase field in order to control the Willmore energies of the level sets individually 
and exclude transversal crossings (which phase fields for De Giorgi’s functional can develop). These 
regularisations may prevent loss of connectedness along a gradient flow in practice, but do not lead 
to a variational statement via T-convergence. 

Furthermore, we would like to emphasise that we can easily describe a weakly* continuous 
evolution of varifolds along which connectedness is lost. Except at one singular time, the varifolds 
are embedded C 2 -manifolds and the evolution is C 2 -smooth. Details can be found in a forthcoming 
paper |Wojl6| . 

It thus is not clear whether the approach of lBil97l does prevent topological transitions, in 
particular, the loss of connectedness, in three ambient space dimensions. At least, it is more 
difficult to implement due to the highly non-linear term including the Willmore energies of level 
sets. 

Also in |Wo j 16| , we show that topological genus is not continuous under varifold convergence 
and that minimising sequences of a constrained minimisation problem with fixed genus may change 
topological type in the limit. 

At this point, we thus know of no other model which can control the topology of phase field 
limits. Furthermore, our results are optimal since they allow us to control as much of the topology 
as can be controlled even for a sharp interface and they allow for efficient implementation. 


1.2. Biological Membranes. Willmore’s and Helfrich’s energies are widely used in the modelling 
of thin elastic structures. The first mention in that context goes back to Sophie Germain |Ger21| . 
Later, their importance has been suggested heuristically |Hel73] based on the principle that when 
a membrane is written as a graph over its tangent space, only derivatives of at most second degree 
should occur in at most quadratic expressions at the base point. Another biological motivation in 
the context of red blood cells is given in [Can70| . 

A Helfrich type functional has also been obtained as a macroscopic limit of certain mesoscale 
models for lipid bilayers 1 PR09 1 LPR14| . Here bilayers are modelled using functions u,v £ 
BV(R n , {0,1}) to express locations of hydrophilic heads and hydrophobic tails of lipid molecules. 
These functions are coupled through a Monge-Kantorovich distance to be close together, con¬ 
strained to satisfy uv = 0 and the perimeter of {it = 1} is penalised. This energy prefers a bilayer 
structure and suitably rescaled versions T-converge to a Helfrich functional with H 0 = 0, xh = 1/2 
and xk = 1/3. 
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A heuristic way to motivate the occurrence of Willmore’s energy in this context comes from 
thin shell theory. In lF.TM02bl lF.TM02al| Friesecke, James and Muller proved T-convergence of 
non-linear three-dimensional elasticity to geometric bending energies including those of Willmore- 
or Helfrich-type in the vanishing thickness limit of thin plates. The admissible class here are 
isometric embeddings of domains in R 2 . The restriction to isometries stems from the fact that 
in-plane stretching energy scales with thickness of the plate and dominates out of-plane bending, 
which scales with the third power of the thickness parameter. Thus in the (second order) vanishing 
thickness limit, we are led to minimise Willmore’s energy in a class of isometric immersions. The 
case of shells (i.e. non-flat structures) has been treated in |FJMM03| . 

Lipid bilayers differ from shells in that they are liquid not solid and as such do not have a 
reference configuration. Assuming inextensibility, we must therefore also minimise over the space 
of Riemannian metrics on the bilayer with fixed area and the isometry constraint turns into the 
fixed area constraint. 

1.3. Further Context. Without any claim of completeness, let us give a bit more context of 
our topic. Willmore’s energy is named for T.J. Willmore who studied it in a series of publica¬ 
tions |Wil651 IWil711 IWil921 IWilOO] and popularised it in his textbook |Wil93| . Independently of 
Willmore’s work, the energy had already been considered as a bending energy for thins plates in 
|Ger21] and as a conformal invariant of surfaces embedded in R 3 in [BT291 ITho23| . 

As mentioned above, Willmore’s energy is famously conformally invariant, in particular scale 
invariant. This shows that both the embedding constraint into and the area constraint 7~L 2 = S 
are needed to give a non-trivial constraint. 

Extrema (in particular, local minimisers) of the Willmore functional are of interest in mod¬ 
els for biological membranes, but they also arise naturally in pure differential geometry as the 
stereographic projections of compact minimal surfaces in S 3 , see e.g. |PS87j . also for examples. 

Even for Willmore’s energy rigorous results are hard to obtain. From the point of view of 
the calculus of variations, a natural approach to energies as the ones above is via varifolds [A1172 ;, 
[HutSbj , where existence of minimisers for certain curvature functionals is proved. The existence of 
smooth minimising Tori was proved by Simon [Sim93j , and later generalised to surfaces of arbitrary 
genus in 111x03. . 

The long-standing Willmore conjecture that W(T) > 47r 2 for all Tori embedded in R 3 was 
recently established in |MN14I , and the large limit genus of the minimal Willmore energy for closed 
orientable surfaces in R 3 has been investigated in [KLSlOj . The existence of smooth minimising 
surfaces under isoperimetric constraints has been established in (Schl2) . A good account of the 
Willmore functional in this context can be found in |KS12j . 

The case of surfaces constrained to the unit ball was studied in [ MR14j and a scaling law 
for the Willmore energy was found in the regimes of surface area just exceeding 47 t and the 
large area limit. While the above papers adopt an external approach in the language of varifold 
geometry, a parametrised approach has been developed in |Rivl4| and related papers. In [KMR14] 
this framework is used to solve the Willmore minimisation problem with prescribed genus and 
prescribed isoperimetric type. I Kill 15 gives a study of the Willmore functional on C 2 -graphs 
and its id-lower semi-continuous envelope. 

Other avenues of research consider Willmore surfaces in more general ambient spaces [LMS111 
ILMIflllMRTl . 

In the class of closed surfaces, if Xk is constant, the second term in the Helfrich functional is of 
topological nature due to the Gauss-Bonnet theorem. So, if the minimisation problem is considered 
only among surfaces of prescribed topological type it can be neglected. The spontaneous curvature 
is realistically expected to be non-zero and can have tremendous influence. It should be noted 
that the full Helfrich energy depends also on the orientation of a surfaces for Hq 0 and not only 
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on its induced (unoriented) varifold. Grofie-Brauckmann |GB93| gives an example of (obviously 
non-compact) surfaces of constant mean curvature H = 1 converging to a doubly covered 
plane. This demonstrates that, unlike the Willmore energy, the Helfricli energy need not be lower 
semi-continuous under varifold convergence for certain parameters. 

Recently, existence of minimisers for certain Helfrich-type energies among axially symmetric 
surfaces under an isoperimetric constraint was proved by Choksi and Veneroni |CV13| . Lower 
semi-continuity for the Helfrich functional on C 2 -boundaries with respect to the ZA-topology of 
the enclosed sets was established by means of Gauss graphs in jPe!97l . 

Results for the gradient flow of the Willmore functional are still few. Short time existence for 
sufficiently smooth initial data has been shown in [SimOll [KS02| (see also [MS03] f and long time 
existence for small initial energy and convergence to a round sphere has been demonstrated in 
[KSOlt lKS12| . Kuwert and Schatzle’s lower bound on existence times in terms of initial curvature 
concentration in space has been generalised to Willmore flow in Riemannian manifolds of bounded 
geometry in |Linl3j . It has been shown that Willmore flow can drive smooth initial surfaces to 
self-intersections in finite time in [| MS03J . This issue seems to be prevented by our connectedness 
functional on the phase field level, although we do not have a rigorous statement on this. Numerical 
simulations suggest that singularities can occur in Willmore flow in finite time |MS02l . i PW071 
Figure 2] gives a numerical example of a disc pinching off to a torus. A level set approach to 
Willmore flow is discussed in [PR04j . 

Studies of numerical implementations of Willmore flow are, for example, due to Garcke, Pziuk, 
Elliott et al. in [BGN081 |Pzi08l 1)11)7 . Particularly interesting here is also an implementation of 
a two-step time-discretisation algorithm due to Rumpf and Balzani [BR12] , 

Phase field approximations of the functional in (1.1), on the other hand, often provide a more 
convenient approach to gradient flows or minimisation of the Willmore or Helfrich functionals. In 
particular if a coupling of the surface to a bulk term is desired, phase field models can provide an 
excellent alternative to a parametric discretisation iPuTOl . 

The idea for applying a phase field approach to Willmore’s energy goes back to Pe Giorgi [ IPG91 !. 
For a slight modification of Pe Giorgi’s functional, reading 


£ 


dG 


(u) = — [ - (s Au —-W'(u)] da; + — f -\Vu\ 2 + - W (u)dx =: W e (u) + X S e (u) 
Co Jn e V £ / Co J n 2 e 


F-convergence to the sum of Willmore’s energy and the A-fold perimeter functional (A > 0) was 
finally proved by Roger and Schatzle (RS06| in n = 2,3 dimensions, providing the lower energy 
bound in the varifold context, after Bellettini and Paolini had provided the recovery sequence 
over a decade earlier [BP93| . First analytic evidence had previously been presented in the form 
of asymptotic expansions in [PLRW05] , This and other phase field approximations of Willmore’s 
energy and their L 2 -gradient flows are reviewed in [BMQ13| . A convergence result for a diffuse 
approximation of certain more general Helfrich-type functionals has also been derived by Belletini 
and Mugnai 1BM101 . 

Regarding implementation of phase field models for Willmore’s and Helfrich’s energy, we refer 


to the work by Misbah et al. in [BKM05] and Pu et al. |PLRW05l IPLW051 IPLW061 IPW071 
PLRW07, PLRW 091 IPulOl IWP07] . In FR W13] . the two-step algorithm for surface evolution has 
been extended to phase fields. A numerical implementation of the Helfrich functional can be found, 
for example, in the work by Campelo and Hernandez-Machado [CHM06j . 


2. Phase Fields 

2.1. General Background. A phase field approach is a method which lifts problems of (n — 1)- 
dimensional manifolds which are the boundaries of sets to problems of scalar fields on n-dimensional 
space. Namely, instead of looking at dE, we study smooth approximations of xe — Xe c - Such 
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models are generally based on a competition between a multi-well functional penalising deviation 
of the phase field function from the minima (usually ±1 instead of the characteristic function values 
0, 1) and a gradient-penalising term preventing an overly sharp transition. Fix f1 <e R ra , open. Our 
model space is 

X = {u £ (H^ n ) | u = — 1 outside f2}= —1 + Wq’ 2 (Q.). 

If dfl £ G 2 , this can be identified with X a = {u £ W 2,2 (VL) \ u = d v u = 0 on <9S7}, otherwise the 
first formulation is technically easier to work with. The boundary conditions are chosen to model 
sets which are totally contained in Q and whose boundaries may only touch dfl tangentially. Then 
(as well as in more general settings) it is known [Mod87j that the Modica-Mortola or Cahn-Hilliard 
functionals 


S E : L 1 (fl) -A- M, S e (u) 


±Ll\Vu\ 2 +lW(u)dx 

+00 


uex 

else, 


with W(u) = (it 2 — l) 2 /4 and Cq = f\ \J2 W(u) d u = 2\/2"/3, T-converge to 


So(u) = ^\Du\(fl) 

with respect to strong .^-convergence for functions u £ i?ld(R n , { — 1,1}) (and +oo else). The 
limit Sq(u) agrees with the perimeter functional Per({u = 1}). The double-well potential W forces 
sequences of bounded energy S e (u e ) to converge to a function of the above form as e — > 0. By 
Young’s inequality |a 2 + eb 2 > ab we can control the i3Y-norms of G{u e ) = /“j \j2W(s) ds 
through 

II G(u £ )||i < C (^J w(u e ) + iy , 

( |VG(u e )|dx = [ \G'(u e )\ \Xu e \dx 
Jn Jq 

< f ^-\G'(u E )\ 2 +e\Xu E \ 2 dx 

Jn 4e 

= 2 S E (u E ). 


So G(u e ) —► v converges strongly in L p (fl) for all p < n/(n — 1) for some v £ BV(Q) by the 
compactness theorem for .BY-functions (up to a subsequence), so in particular pointwise almost 
everywhere (further subsequence). ||ue|||, 4 (f 2 ) < C (l + j f n W{u E )dx), so u E —> u weakly in L 4 (f2) 
and u E —» G _1 (u) pointwise a.e. which together implies that u E -A u = G“ 1 (u) strongly in L p (fl) 
for all p < 4. The functional has originally been studied in the context of minimal surfaces, see 
e.g. [Mod871 ICC951 ICCOfi) . We will view them in a slightly different light. If u E is a finite energy 
sequence, we denote the Radon measure given by the approximations of the area functional as 

He{U) = - f A |Vu e | 2 + - W(u E ) dx. 

Co Ju 2 e 

Taking the first variation of this measure and integrating by parts, we get 

<5/+(<(>) = — [ (-e A u e + - W'(u e )) cj)dx. 
c o Jn V e J 


The integrand 

v e := -e A u e + - W'{u e ) 

£ 

is a diffuse analogue of the mean curvature as the variation of surface area. A modified version 
of a conjecture by De Giorgi then proposes that squaring v E , integrating over SI and dividing by 
e should be a T-convergent approximation of W(dE) for G 2 -regular sets E, again with respect to 
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strong /^-convergence. Dividing by e is needed to convert the volume integral into a diffuse surface 
integral over the transition layer with width of order e. This is true at least in n = 2,3 dimensions 
|RS06I . From now on, we will restrict ourselves to these dimensions. We thus introduce the diffuse 
Willmore functional 

, . f y- Jo y (eAu - i W'{u)) 2 dx u € X 

WJu) = < c ° JQ e v £ 

|+oo else 

and the Radon measures associated with a sequence it e with bounded energy S E (u E ) + W E (u E ) 

Vs(U) = — ( '-\S7u e \ 2 + -W{u E )dx ) a e {U) = — f-v 2 dx. 

Co Ju 2 £ Co _/[/ £ 

Since the measures /z e ,a e are uniformly bounded, we can invoke the compactness theorem for 
Radon measures to see that there are finite Radon measures p, a with support in Q such that, up 
to a subsequence, 

u E —» u strongly in L 1 (f 2 ), p E /z, a E A a _ 

A finer comparison shows that \Du\ < 2p. Clearly by construction u = Xe — Xe° for some E C f2, 
and from the above argument is not difficult to see that E is a Caccioppoli set. In three dimensions, 
however, E may well be empty. To see this, one can use the construction of |MR14j . where two 
concentric spheres of very similar radius are connected by a modified catenoid. As the catenoid is 
a minimal surface, this keeps the Willmore energy bounded. Taking the radii to converge to the 
same quantity and approximating them with the recovery sequence of a phase field, we see that E 
is the L 1 -limit of the enclosed set between the two original sphere, so empty. 

In [RS06] it has been shown that p is in fact the mass measure of an integral n — 1-varifold with 
| d*E\ < /i and H 2 ■ fj, < a for the generalised mean curvature H^ of fi. Unlike /z, a can still behave 
wildly. 

We aim to minimise the Willmore energy among connected surfaces with given area S. To 
prescribe area /x(U) = S , we include a penalisation term of 

(S. - S) 2 

in the energy functional. So we have a good approximation of Willmore’s energy available, the 
area constraint can be approximated through an easy penalisation, and confinement to Q comes 
automatically from the choice of function space X. The main challenge left is controlling the 
connectedness of spt(/z) which we interpret as a sharp interface membrane. 

2.2. Equipartition of Energy. Intuitively, the Cahn-Hilliard energy prefers phase fields to be 
close to ±1 in phase and make transitions between the phases along layers of width proportional to 
e. In particular, to attain the minimum, equality |Vit e | = IG^Ue)! is needed in Young’s inequality 
on most of ft. So it makes sense that the terms | |Vu E | 2 and 7 W(u e ) should make an equal con¬ 
tribution to the measure p, £ . The difference of their contributions is controlled by the discrepancy 
measures 

UU) = ~ [ x |Vu e | 2 — -W(u E )dx, 

Co Ju 2 £ 

their positive parts £ £i+ (t/) = A (| |Vu E | 2 — - W(u e )) + dx and their total variation measures 
|£ e |(U) = A || |Vrt e | 2 — ^W(u e ) | d.T. We will see that for a suitable recovery sequence, |£ £ | 
vanishes exponentially in e; more generally the discrepancy goes to zero in general sequences u E 
along which a e + /+ stays bounded as shown in l iSOli . Propositions 4.4, 4.9]. The result has 
been improved along subsequences in [BM10 , Theorem 4.6] to L p -convergence for p < 3/2 and 
convergence of Radon measures for the gradients of the discrepancy densities for an approximation 
of Helfrich’s energy. 
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The control of these measures is extremely important since they appear in the derivative of the 
diffuse (n — l)-densities r 1 ~ n p e (B r (x)) with respect to the radius r. 


2.3. Connectedness. In order to ensure that the support of the limiting measure p is connected 
we include an auxiliary term C e in the energy functional. The heuristic idea behind this is that if 
the support of the limiting measure spt(/x) is connected, then so should the set {pi < u e < P 2 } for 
— 1 < pi < P 2 < 1. These level sets away from ±1 can be heuristically viewed as approximations 
of spt(p), and in other situations, they can be seen to Hausdorff-converge to it |CC95j . This is not 
quite true in our situation (in three dimensions, see Remark |2.8[ ), but the intuition still holds. 

Our concept is to introduce a quantitative notion of path-connectedness and penalise the mea¬ 
sured disconnectedness. Take a weight function F £ C°(M) such that 

F> 0, F = 0 on [pi,p 2 ], F(-l), F(l) > 0. 


and the associated geodesic distance 


d F ^ u ' ) (x, y) = inf |y F(u) dH 1 


K connected, x,y € K : 'H l {K) < u{e) >, w(e) —> 00 as £ —> 0 


In particular, w(e) = 00 is not excluded. If {p\ < u e < P 2 } is connected, we can connect any two 
points x, y £ u“ 1 (pi, p 2 ) by a curve of length zero. If it is not, then d F ^ u \x, y) gives a quantitative 
notion of how badly path-connectedness fails between these two points. To obtain a global notion, 
we take a second weight function <f> £ C c (— 1,1) resembling a bump, i.e., 

4> > 0, {</> > 0} = (pi,p 2 ) (-1,1), J 4>(u)du>0 


and take the double integral 

C E (u) = \[ [ <t>{u{x))^{u{y))d F{u) {x,y)dxdy. 

£ Jn Jn 

So, if {<f>(u e ) > 0} = {p\ < u e < P 2 } is connected, we can connect any two points x,y £ such that 
<j)(u e (x)) <t>{u e (y)) > 0 with a curve of length zero, hence d F ( Ue \x,y) = 0 and both the integrand 
and the double integral vanish. 

If on the other hand spt(/z) is disconnected, then we expect that d F ( Ue ' ) should be able to discern 
different connected components such that liminf £ _,.o C e {u E ) > 0. The core part of our proof is 
concerned with precisely that. We need to show that <j> detects components of the interface and 
that d F ^ Ue ^ distinguishes them. For the first result, we need to understand the structure of the 
interfaces converging to p, and make sure they cannot be so steep in u~ 1 {p\ 1 P 2 ) that the double 
integral does not see them in the limit. 

For the second part, the challenge is to understand how phase fields converge away from the 
interface. It is clear that they converge pointwise almost everywhere, but it is not a priori clear 
whether they can stay away from ±1 along a curve in three dimensions. Namely, if the support of 
a limiting measure p of p E had two components C \, C 2 , but the approximating sequence satisfied 
u e £ (pi, P 2 ) in a neighbourhood of a curve connecting C\ and C2, then d F ^ Ue ^ could not 
distinguish between the components. 

In two ambient dimensions, this is excluded relatively easily. In three-dimensional space however, 
the curve would not be seen by the limiting measure p , which is an integral 2-varifold. Thus 
we need to show that in analogy to the sharp interface, a long and narrow tunnel-like strucure 
around blows up in diffuse Willmore energy. We like to think of this as a convergence localised 
on curves, i.e. sets of co-dimension n — 1 = 2. 
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2.4. Main Results. We have the following results for phase fields with bounded Willmore energy 
and perimeter. In the following we assume that u E denotes the continuous representative. 


Theorem 2.1. Let n = 2 and u e £ X a sequence such that 

limsup (We + S e ) (u e ) < oo. 

£—>-0 

Denote p = lim e _>.o Me for a subsequence along which p E converge weakly as Radon measures to a 
limit and take any Q,' <e R n \ spt(p). Then |it e | —> 1 uniformly on 12'. 

This result has already been proved by Nagase and Tonegawa in jNT07j . where they also show 
a r — lim inf-inequality for W E + S e in two space dimensions. The use of a monotonicity formula in 
their proof resembles ours with differences in the estimate on the discrepancy measures. We show 
our proof here to illustrate some ideas in the more involved argument for n = 3. 


Theorem 2.2. Let n = 2 and D <g (—1,1). Then u e 1 (D) —> spt (p) in the Hausdorff distance. If 
n = 3, up to a subsequence, uf 1 (D) Hausdorff-converges to a set K which contains spt(/i)- 


In three dimensions, K ^ spt(/i) is possible, see Remark 2.8 As an illustration, the following 


statement describes the convergence of phase fields on lower dimensional objects that we need. 
In the proofs a slightly different version of the theorem is used to account for changing curves as 
£ -A 0. 


Theorem 2.3. Let n = 2 and u e £ X a sequence such that 

limsup (}V e + S e ) (u e ) < oo. 

£—>-0 

Denote p = lim^o p E for a subsequence along which p E converges weakly as Radon measures to a 
limit. Take any K <s R n \spt(M) such that 0 < 'H 1 {K) < oo. Then \u E \ —> 1 converges TL 1 \K~almost 
everywhere. 

For our application to connectedness, we define the total energy of an £-phase field as 


, 9n W E (u) + e-°(S £ (u)^S) 2 +e-«C E (u) u£X 

(2.1) £ E (u) = < 

I +oo else 

for a £ (0,4), k > 0. 

Remark 2.4. If w(e) < oo, existence of mininimisers for the functional £ e is a simple exercise in the 
direct method of the calculus of variations, since uniform convergence of a minimising sequence 
u e ,k £ W 2 ’ 2 (f2) -A C°({2) for fixed £ guarantees convergence of the distance term. 


Theorem 2.5. Let u e £ X be a sequence such that (W e + S E )(u E ) < C for some C > 0 and p, a 
Radon measures such that p E p, a e a. If spt (p) is disconnected, then 

lim inf C e (u E ) > 0. 

£—^0 

Using these results, we can show the following. 

Theorem 2.6. Let n = 2,3 and u E £ X a sequence such that liminL-^o£ e(w e ) < oo. Then the 
diffuse mass measures p E converge weakly* to a measure p with connected support spt (p) C U and 
area p(Ll) = S. 

Using IlSOtil . p is also the mass measure of an integral varifold. The main result of I RS06] can 
be applied to deduce T-convergence of our functionals in the following sense: 
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Corollary 2.7. Let n = 2,3, S > 0, Q <e R ra and E <g LI, with smooth boundary dE £ C 2 with 
area TL n ~ 1 ( y dE) = S. Then 


rC^Hfi)) - lim £c(xe - Xe°) 


W(dE) dE is connected 
+oo otherwise 


Here we have adopted the notation of |R.S06| where T-convergence is said to hold at a point 
if the lim inf- and lim sup-inequalities hold at that point. This distinction is necessary since it is 
not clear what T-convergence properties hold at other points u £ BV(Sl, { — 1,1}) since W E does 
not converge to the If -lower semi-continuous envelope of W at non-embedded sets even in two 
dimensions. 

This issue stems from the fact that the stationary Allen-Cahn equation —A u + W'(u) = 0 
admits global saddle-solutions which have vanishing sets asymptotic to a union of lines. The 
best known example vanishes along the coordinate axes and is positive in the first and third 
quadrants and negative in the second and fourth ones, see jdPKPWlOl . These solutions can 
be used to approximate transversal crossings with zero Willmore energy, while the lower semi- 
continuous envelope of Willmore’s energy/Euler’s elastica energy becomes infinite at those points. 
Interestingly enough, the crossing has zero Willmore-energy as a varifold, so that the energy can 
still be justified, despite the fact that it does not give the sensible result in our situation. 

While heuristic considerations suggest - and numerical simulations appear to confirm - that at 
least in n = 2 dimensions and when 0 ^ spt (</), our additional term in the energy might prevent 
saddles, we do not investigate convergence at non-embedded points further. 


Remark 2.8. Uniform convergence generally fails in three dimensions. This can be seen by taking 
a sequence u e as the optimal interface approximation which will be given in the proof of Corollary 
2.7 and adding to it a perturbation g((x — Xq)/e) for any g £ C“(R n ). It can easily be verified 


that the energy will remain finite (at least if \g\ -C 1) but clearly the functions do not converge 
uniformly. 

Heuristic arguments show that for modifications u E = u e +e“ g{e~ !i ■), the coefficients a = /3 = 1 
should be optimal in three dimensions and a = 1/2, /3 = 1 if n = 2, so we expect a convergence 
rate of yfe in the two-dimensional case. 

Since uniform convergence does not hold, it follows that the pre-images of certain level sets can 
have additional points in the Hausdorff limit. 


Remark 2.9. Since S e (u e ) is uniformly bounded, a standard argument for the Modica-Mortola 
model shows that u e -A \e — Xn\E i n L p for p < 4/3. We are going to show in Lemma 3.1 that 
the sequence u e is uniformly bounded in L°°(W l ) when additionally limsup^g W(u e ) < oo, which 
implies L p -convergence for all p < oo. 


3. Proofs 


The proofs are organised in the following way. First, anticipating the results of section 3.4 
we show T-convergence of £ e to W. We decided to move the proof to the beginning since it is 
virtually independent of all the other proofs in this article and can be isolated, but introduces the 
optimal interface transitions which will be needed later. After that, we proceed chronologically 


with technical Lemmata (sections 3.2 and 3.3) and the proof of Theorems 2.1 2.3 
in section [+T1 


2.2 2.5 and 2.6 


3.1. Proof of 1 -Convergence. We now proceed to prove Corollary |2.7| 


Proof of the liminf -inequality: It will follow from Theorem 2.6 that £ e (u e ) —> oo if dE is discon¬ 
nected. If dE is connected, the main part of this inequality is to show that if u e —> xe — Xe° 
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in L 1 (fi) and <5 + 1, then liminfe^o £e(u E ) > W(dE). Since £ e > W e and enforces the 

surface area estimate, we obtain with 1RS06I that 


liminf £ e (u e ) > W(dE). 


□ 


Proof of the limsup -inequality: We may restrict our analysis to the case of connected boundaries 
with area T~L n ~ 1 i y dE ) = 5. In this proof, we will construct a sequence u E £ X such that 

lim ||u £ - (xe - XE=)lli,n = 0, lim£ e (u £ ) = W(dE). 

e-10 e->0 

The construction is standard and holds for arbitrary n £ N. We take the solution of the one¬ 
dimensional problem and apply it to the signed distance function of dE , thus recreating the shape 
of dE with an (approximate) optimal profile for the transition from —1 to 1. For <5 > 0 consider 

XJs '■= {x £ R n | dist(a;,9.E) < <5}. 

Since E (e ft, Us C fl for all sufficiently small S, and since dE £ C 2 is embedded, there is 5 > 0 
such that 

if : dE x (—5 , 6) -» Us, ip(x, t) = x + tu x 

is a diffeomorphism. We have the closest point projection 7 r : Us —> dE and the signed distance 
function d(x) = sdist(x, dE) from dE which is positive inside E and negative outside E. Both are 
C 2 -smooth on Us and they satisfy if(n(x), sdist(x)) = x and (see e.g. jGTOll Section 14.6]) 

Vd(x) (x); Z\d(x) ' d(x) + 0(d(x) ) 

where C n ( x ) is a constant depending on the base projection point n(x). Let us consider the optimal 
transition between —1 and 1 in one dimension. This optimal profile is a stationary point of 5i, i.e. 
a solution of —q" + W'(q) = 0 (and thus a zero energy point of Wi) with the side conditions that 
linit^ioo q(t) = ±1. Note that the optimal profile satisfies 

q" = W\q) =► q"cf = W'{q)q' =► ^ {{q'? - W(q)) = 0, 

so we find that ( q') 2 — W(q) = c. We look for transitions where both {q') 2 and W(q) are integrable 
over the whole real line, and since limt-^oo W(q(t)) = 0, we see that c = 0 and (q') 2 = W(q). This 
gives us equipartition of energy already before integration. For simplicity, we focus on W(q) = 
( q 2 — l) 2 /4 which has the optimal interface q(t) = tanh(t/-\/2). Note that the functional rescales 
appropriately under dilations of the parameter space so that we have equipartition of energy before 
integration also in the e-problem. The disadvantage of the hyperbolic tangent is that it makes the 
transition between the roots of W only in infinite space, so we choose to work with approximations 
q e £ C 00 ^) such that 

(1) q E (t) = q(t) for |t| < 8/{3e), 

( 2 ) q e (t) = 1 for t > <5/( 2e), 

(3) q e (-t) = —q s (t), 

(4) q’ E > 0 , 

(5) \q”\(t) < e 3 for all |f| > S/e. 

Such a function is easily found as tanh converges to ±1 exponentially fast. Then we set 

u £ (x) = q e {d{x)/e). 

A direct calculation establishes that with this choice of u e , we have |5 e (n e ) — 5| < e 7 for all 7 < 2, 
ICel < £m f° r all m £ N and lim e ^. 0 W e (u e ) = W(dE). It remains to show that lim^o e~ K C e (u E ) = 
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0. We will show that even C E (u e ) = 0 along this sequence. Since dE is connected and ip is a 
diffeomorphism, all the level sets 

{u e = P} =i’(dE,eq~ 1 (p)) 

are connected manifolds for p £ (—1,1) and e small enough such that q E = q around q^ 1 (p). We 
know that 

{(p{u e ) > 0} = {pi <u e < p 2 ) 

and pick any p £ (p\,p 2 ). Now let x,y £ fl, cp(u £ (x)), cp(u e (y)) > 0. We can construct a curve 
from x to y by setting piecewise 


7i : M(a:)]] ->• ^ , 7i {t) = n(x) + tv 
73 : [0, (%)] -A n , 73(i) = n(y) + t u^ y) 

and 72 any curve connecting tt(x) to ir(y) in {u e = p}. This curve exists since connected manifolds 
are path-connected. The curve 7 = 73 © 72 © 7] -1 connects x and y and satisfies by construction 
<^>(7 (t)) > 0, so F( 7(i)) = 0. Therefore we deduce 

d F ^\x,y) = 0 


if <p{u e (x)), <p(u e (y)) ^ 0 since the connecting curves have uniformly bounded length and ui(e) —> 00 . 
Thus in particular 

4 [ (p{u e (x))(/)(u e (y))d F{ue Xx,y)dxdy = 0. 

£ Jnxn 


□ 


3.2. Auxiliary Results I. A lot of our proofs will be inspired by 1RS06] which again draws from 
[HT00] . We will generally cite [RS06J because it treats the more relevant case for our study. In 
this section, we will prove a number of auxiliary results which concern either general properties 
of phase fields or properties away from the support of the limiting measure p, which will enable 
us to investigate their convergence later. Results concerning phase interfaces are postponed until 
section 13.31 

We start with an optimal regularity lemma for phase fields on small balls. This is an improve¬ 
ment upon [RS061 Proposition 3.6] where in three dimensions, only boundedness in L p for finite 
p and a slow growth in L°° could be obtained. The new i°°-bound also implies local Holder 
continuity which is used extensively later. 

Lemma 3.1. Let n = 2,3, H <g M" and u e £ X such that a := limsup e _ >0 W e (u £ ) < 00 . Then 
there exists £q > 0 such that for all e < £q 

(1) ||'u £ ||oo,R ’ 1 — C where C is a constant depending only on a and n. 

(2) u e is 1/2-Hdlder continuous on £-balls, i.e. 

C 1 

I U e {x) - u E (y )I < —j= \x - yX V x £ K", y £ B(x,£). 

V£ 

Again, the constant C depends only on a and n. 

Proof. We will argue using Sobolev embeddings for blow ups of u E onto the natural length scale. 
In the first step, we show the set {|u 6 | > 1} to be small. In the second step, we estimate the 
L 2 - norm of the blow ups, in the third we estimate the full W 2 ’ 2 -norm and use suitable embedding 
theorems to conclude the proof of regularity. 
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Step 1. First observe that 

a £ ( 0 ) > a £ ({u e > 1 }) 

1 

c o J{u e > 1} 


J j ^ W'(u e )^ dx 

= —[ w\u e )d vUe <m n - 1 

CO £ Jd{u £ > 1} 

+ — f e(Au £ ) 2 + -W"(u £ )\Vu £ \ 2 + \(W'(u e )) 2 dx 

C-0 J {Ur >1\ ^ ^ 


( 3 . 1 ) 


1 


> — f e(Au £ ) 2 + 4 |Vu £ | 2 + \ {W'(u £ ) ) 2 dx 

CO J {u e >l} C £ 

using that W”(t) > 2 for t > 1. The boundary integral vanishes if {u £ > 1} is of finite perimeter 
since u £ € W 2 ' 2 c —► C 0,1 / 2 is continuous and W'(l) = 0. If this is not the case, take 9 \ 1 
converging from above such that {u £ > 9 } is of finite perimeter. This holds for almost all 9 £ R. 
The sign of the boundary integral can be determined since d„u £ < 0 on the boundary of {it £ > 9} 
and W'(9) > 0 for 9 > 1 so that the same inequality can still be established. By symmetry, the 
same argument works for {u £ < —1}. 

Step 2. Let x e £ Q be an arbitrary sequence and define the blow up sequence u £ : K" —> K. by 

u £ {y) = u £ {x e +Ey). 

Then we observe that 


[ u 2 e dx= [ ( Itt e I - 1 + l ) 2 d y 

JB ( 0 , 2 ) JB( 0 , 2 ) 

< [ ( (\u £ \ - 1 )+ + l) 2 d y 

JB( 0 , 2 ) 

<2 [ (|w £ | - 1)| + 1 dy 

Jb( 0,2) 


< 2 e 3 ~ n / -3 W'(u £ ) 2 dy + 2 " +1 u n 

{\ u e I >1} £ 

<2 {2 n u n + c 0 E 3 - n a £ (ti)} . 

As usual, oj n denotes the volume of the n-dimensional unit ball. In exactly the same way with a 
slightly simpler argument we obtain 


' B(0,'2) 


(W'(u e )) 2 dy<C(a,n). 


Step 3. Now a direct calculation shows that 

[ (A u e - W'(u £ )f dy=[ (e 2 Au £ - W’(u E )) 2 ( x e + ey) dy 

JB( 0,21 J B( 0 , 2 ) 


>B( 0 . 2 ) 

= c 0 E 3 ~ n a E {B(x e , 2 e)). 


Thus 


A U e 11 2 , T? (0,2) < || A U E — W\u £ ) 11 2 , A? (0,2) + 11 W'(u E ) ||2,B(0,2) 
< \Jc 0 s 3 ~ n a e (Q.) + y/C(a,n). 


I u £ 1)2,15(0,2) + II A u e H2,s(o,2) < C(a, n) 


In total, we see that 
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for all 0 < e < 1 so small that a s (Q) < a(Q) + 1. Using the elliptic estimate from 1(1 1011 Theorem 
9.11], we see that 

ll^e|| 2 , 2 ,B( 0 ,l) < C{a,n), 

where we absorb the constant depending only on n and the radii into the big constant. Using the 
Sobolev embeddings 

W 2 ' 2 {B( 0,1)) W ll 6 (B(0,1)) 


we deduce that 

|we|o,i/ 2 ,B(o,i) ^ C(n,a), 

again absorbing the embedding constants into the constant. In particular, this shows that 

||«e||oo,B(o,i) < C(n,a). 

But since this holds for all sequences x E , we can deduce that 

11 11 e | |oo,M™ fs C{jl, Cxfj. 

Furthermore, for x, z G R n with \x — z\ < e < £o, we choose x e = x to deduce 

\u E (x) - u E (y) | = |u e (0) - u E ((y-x)/e)\ < C(n, a) \ (y-x)/e | 5 = \ x - y \h, 

v e 


□ 


Remark 3.2. Without prescribing boundary conditions as in our modified space X, the result 
could still be salvaged on compactly contained subsets. Techniques for estimating quantities over 
{Kl > 1} in that case can be found in (RS06i Proposition 3.5], which we include below for the 
readers’ convenience. 

Proposition 3.3. iTSTili. Proposition 3.5] For n = 2,3, fl C R n , e > 0, u E G C 2 (f2), v E G C°(U), 

—eAu e H— W'(u E ) = v e in fl, 

£ 

and SI' g !1, 0 < r < dist(U', dt2), we have 

f W\u E ) 2 < C k ( 1 + r~ 2k e 2k )e 2 f v 2 + C k r~ 2k e 2k [ W'(u e ) 2 

J{|u e |>i}nn' In J{|u c |>i}nn 

for all k G No- 

Thus u e is uniformly finite on compactly contained sets for small £ and steps 2 and 3 from the 
proof of Lemma [rT] go through as before. A useful rescaling property is the following observation 
from the proof of ]RS06| Theorem 5.1]. 

Lemma 3.4. Let u e : B(x,r) —> R, A > 0 and u E : B(0,r/X) —► R with 

u E (y ) = u e (x + X y). 


Set f := r/X, e := e/X, 

£ e :=- ( £ -\S7u e \ 2 + \w{u E )]c n , a E :=—feAu E -\w , (u E ))c n . 

Co \2 e J Co s \ e ) 

Then 

r 1 ~ n fi e (B( 0, f)) = r 1 ~ n p e {B(x, r)), f 3 ” n a e {B( 0, f)) = r 3 ~ n a e {B(x, r)). 

The discrepancy measures behave like p under rescaling. 

Proof. This can be seen by a simple calculation similar to the one in the proof of Lemma |3.1| □ 

For the reader’s convenience we include the following classical monotonicity result. 
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Lemma 3.5. [ RS06t Lemma 4.2] For x £ K" we have 


dp ^ y p n c 0 p n+1 


1 


Co P JB(x,p) 


' dB(x,p ) 

v E (y - x , Vu £ ) dy. 


e{y - x^Ue) 2 dU n 1 {y) 


In low dimensions n = 2,3, the second and third term in the monotonicity formula can easily 
be estimated after integration. While the result is known, we fixed some details in the proof of 
I lSOli. Proposition 4.5], so we include it here for completeness. 


Lemma 3.6. [ RS061 Proposition 4.5] Let 0 < r < R < oo if n = 3 and 0 < r < R < 1 if n = 2, 
then 


(3.2) 


r 1 -' n p e (B(x, r)) < 3 R^^B^x, R)) + 2 f &■+(£& P)) dp 

Jr P 

^ j.3 —n 

+ 2{n- f )2 R )) + (—ijj a e (B{x, r)) 

o2 ol-n 


where Rq := minji?, _Rq} and i?n is a radius such that f l C B( 0, i?n/2) 

Proof. Without loss of generality we may assume that x = 0 and write B p := B(0,p), f(p) = 
p 1 ~ n y e (B p ). Observe that for any function g : Br M. we have 




g(x) dec dp 


■>b r 

1 

n— 1 


50 ) 


/ max{|a:|,r} 


p n dp dx 


’b r 




max{|x|, r} 


n—1 


i 

i ?"- 1 


dx 




<?(x) dT-T" 1 


d p= f 

JB R \B r 



dx. 


and 
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Using this to integrate the derivative we obtain using Young’s inequality 

i-R 


f(R)-f(r)= [ >(p)dp 
J r 




e{\7u e ,y) 2 1 v e (y,S7u e ) 


coJ BR \B r \y\ n+1 n-l |p|’ 


zr— dp 


1 


> 


(n- 1 )c 0 r n - 1 J B 

r R -U+(B p ) 


T / v e (y, Vile) d y - 
Jb„ 


(n - 1) c 0 R n 1 J B] 


Ve (p, Vn e )dp 


1 


dp 

e (Vu e , p) 2 1 


co 


+i 


(y,v^ e ) 2 , i 2 \ , 

n - 1 V |y|2(n—1) + 4 („ _ 1) £ V e) 


x s(y,Vu s ) 2 + |p| 2 ug 


cq 7* n—1 J Br 2|p| 2 2X(n — l) 2 £ 


X 


e{y,Vu e ) 2 1 |p| 


+ 


coi?- 1 2 |p| 2 2A(n-l) 2 e 


dp 


> 


£e,+ (-®p) 


dp — 


4 (n — l) 2 JB R \B r £ 


v 2 dy 


-A/(r)- — 


2A (n — l) 2 r n ~ l J Br £ 


— dp 


-A f(R)~ 


R 2 


1 


2A (n — l) 2 i?" _1 J Br £ 


dp. 


where A £ (0,1). Here we used that n = 2,3 to obtain that 2(n — 1) < n + 1, so that |p|" +1 < 
|p| 2 (” _ i) for all \y\ if n = 3 and for |p| < 1 if n = 2. When we bring all the relevant terms to the 
other side, this shows that 


(1 + A) f(R) — (1 — A) f(r) > — [ R dp - 

Jr P 


P" 

3 —n 


4 (n — l) 2 


a e (B R \ B r ) 


r (B r ) - 


R 2 


2\{n-l) 2 " 2A (n - l) 2 i?"- 1 

Setting A = 1/2 and multiplying by two proves the Lemma. 


r (Br). 


□ 


Remark 3.7. If n = 3, we may let R -A oo and subsequently e —> 0, r —> 0 and finally A —> 0 that 
we have 

limsupr 1 ""p(H(a;,r)) < 1 a(fl) 

r-> o 4 (n-l) 2 

at every point x £ M 3 such that a({x}) = 0 (i.e. when lim r _>.o a(B r ) = 0). Using the results of 
[ RS061 . p, is an integral varifold, so this yields a Li-Yau-type ILY82I inequality 

(3.3) 9*{p,x) = limsup l Jj (B( x d r )) <- __j_Q,(^), 

r _>o nr 1 d7t 

This inequality is usually found with a 4 in place of the 16 which stems from a different normali¬ 
sation of the mean curvature and W(p) in the place of a, see also the proof of |Top98[ Lemma 1] 
or |KS12[ Proposition 2.1.1], 
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In n = 2 dimensions, we had to assume R < 1. Indeed, an inequality of this type cannot hold 
since circles with large enough radii have arbitrarily small elastic energy. Still, setting R = 1, a 
similar bound on the multiplicity in terms of a and S can be obtained. 


(3.4) 


The version of the monotonicity formula (3.2) which we will use is the simplified expression 

rR ^,+ (B{x,p)) 


r 1 n p e (B(x,r)) < 3R 1 n p e (B(x, R)) + 3 a e (B(x, R)) + 2 


dp. 


This holds generally if n = 3, and when R < 1 if n = 2. Furthermore, we have the following estimate 
for the positive part of the discrepancy measures. It is a precise quantitative refinement of the 
classic statement that smooth solutions of the stationary Allen-Cahn equation — Am + W'(u) = 0 
on R" satisfy |Vu| 2 < 2 W(u) |Mod85| . 


Lemma 3.8. i RS061 Lemma 3.1] Let n = 2,3. Then there are i5 0 > 0, M G N such that for all 
0 < <5 < <5o > 0 < e < p and 

po := max{2,1 + S~ M e} p 

we have 

P 1 ~ n tc,+ (B(x : p)) <CS p l ~ n p. e (B( x, 2 p)) + C 5~ M £ 2 p 1 "" f - v 2 da; 

JB(x,p 0 ) £ 

+ CS- M £ 2 p 1 ~ n [ \w'{u E ) 2 dx + — . 

aS(x,po)n{|u e |>l} e P 


The following observation is the key ingredient in order to obtain our required convergence 
results. Geometrically, it expresses that phase fields u e which stay away from ±1 at a fixed 
point need to have a minimal positive measure a or p. In three dimensions, this is only proved 
under additional technical assumptions. The proof uses infinitesimal Holder continuity and the 
monotonicity formula to transition from the c-length scale to the macroscopic scale. 

In a slight abuse of notation, we will denote the functionals defined by the same formulas by 
W e , S e again, although they are given on spaces over B\ := B( 0,1) instead of H. 


Lemma 3.9. Let n = 2,3, 9 £ (0,1), 0 < r\ < 1/2. Consider the subsets 

Y 2 :={uelf 2 ’ 2 (Bi) : ]u(0)| <0} 


in n = 2 dimensions and 

Y e 3 := J u G W 2 ' 2 {B i) : |u(0)| < 9 and a e [B+ [ W_^f_ ^ ^ ^ 

[ JB^n{\u s \>i} £ 

Define T e : IF 2,2 (I3 1 ) —> [0, oo) as 


R e {u) = We (it) + S s (u). 

Then 9 q ;= liminfg-^o inf ue y 2 T e (u) >0 if n = 2 and 

9q := liminf inf JF^(u) > 0 

ueY? v ’ 

if n = 3. The same works if instead u(0 ) >1/9. 


Proof: By Holder continuity, the condition that |it(0)| < 9 leads to the creation of an infinitesimal 
diffuse mass density £ 1 ~ n p e (B e ) > c n ^ & ,e- We will use the monotonicity formula to integrate this 
up to show that if a = 0, macroscopic mass is created as well. In three dimensions, there is an 
additional technical complication which forces us to make two steps, one from the e-scale to the 
length scale of y/e and a second one to the original scale. 

Step 1. In a first step we show that the diffuse mass densities are uniformly bounded on large 
enough length scales. For i£ll, set f e {p ) := p 1 ~ n p e (B(x, p)). Without loss of generality, we may 
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assume that f E ( 1) < 1 for small enough e > 0 and that a e (B 3 /^) + / B n { \ u | >x j j:jW'(u E ) 2 dx < 1 
due to [R.S061 Proposition 3.5] ( incl uded here as Propositio n |3.3[ ). 

Take 5 = log(e:) -2 in Lemma 3.8 to obtain from Lemma [3.6| that for 0 < r < R = 1 we have 


fe(r)<3f £ (R) + 3a e (B R ) + 2 ^ + ^ p) dp 

Jr P 


r R 


< Ca.n + c 


1 fe(2p) 


log(e) 


log(e) 2 p 2 


+ e 2 log(e) 2M p n ^ a e (B p ) + j 


s p n{|u £ |>i} 


d. 


C 


r R 


— C a n T . . \y I 

log(e) 2 J r 
C e 2 log(e) 2M 


fe(2p) 


dp - 


Ce 


< Ca.n + 


n — 1 
C 

log(e) 2 


P log(e) 2 

- (r 1_ra — R l ~ n ) 

/e(p) 


1 1 

r i? 


P 


■dp 


for a uniform constant C a , n for r > e if n = 2 and for r > i/e if n = 3. We use Gronwall’s 
inequality backwards in time to deduce that 


/ (j r 2 R i \ 

Mr) < C a , n exp l ~ dpj < C & , r 


on [e, 1) if n = 2 and on [y/e, 1) if n = 3. 
Step 2. If n = 3 and additionally 

MB^) 


/s^n{|«el>i} 




we can estimate the terms in the second line more sharply to obtain uniform boundedness of f e (r) 
also for e < r < R = y/e. 

Step 3. Now we turn to the proof of the statement. For a contradiction, assume that ( a e + 
p E )(B i) —> 0 for a suitable sequence u E . The functions u E are C 0,1//2 -Holder continuo us w ith Holder 
constant C/y/e for a uniform C > 0 on Bi/ 2 , as can be obtained like in Lemma 3.1 The only 
difference is that we need to use FRS06t Proposition 3.5] (included here as Proposition 3.31 to 
estimate /n„ j >:l i ^sW {u e ) 2 dx due to the lack of boundary values. It follows that u{x) < ^ for 
x £ B(0,ce) for some uniform c > 0, which implies 


£ Jb( 0,ce) V 2 J 


dx =: c n 


> 0. 


In the following, we will assume n = 3. The two-dimensional case follows with an easier argument 
of the same type. First we deduce that 

fs( £ ) < 3 f E (y/e) + 3 a e (B^) 

r^ e 1 e - - - 


C 


log (e) 2 P log(e) 2 p : 


+ £ 2 l0g(£) 2M p-"£^dp. 
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Since the second line goes to zero as e —> 0 and a e —*• 0 by assumption, we deduce that f E {y/e) > 
c n ,a,e /4 for all sufficiently small e > 0. Finally, we obtain 


fe(Vz) < 3/ e (l) + 3a e (Bi) 


C 


VE !og {e) 2 P 


log(e) 2 p 1 


Fe 1 log {eY M p- n dp. 


Again, the terms in the second line vanish with e -> 0 and a e (B i) —> 0 by assumption, thus 


fe( 1) > 


Cn,cx,6 

16 


for all sufficiently small e > 0. But this contradicts the assumption that (a £ + p e ){Bi) —> 0, so 
we are done. The case it(0) >1/0 follows similarly in two dimensions; in three dimensions it is 
automatically excluded by the choice of space Y/ 3 . □ 


The estimate above is not sharp, but suffices for our purposes. After this Lemma, everything is 
in place to show uniform convergence of phase fields in two dimensions in section |T4} while further 
results will be needed for Hausdorff convergence of the transition layers. In a very weak phrasing, 
Lemma 3.10 suffices, more precise versions (and the application to connectedness) need the entire 
section mr 


3.3. Auxiliary Results II. In this section, we will derive technical results concerning how phase 
field approximations interact with the function </> as needed for the functional C e to impose con¬ 
nectedness. While the previous section focused on estimates away from the interface, here we 
investigate the structure of transition layers close to spt(/i). The following Lemma is a special case 
of |RS06t Proposition 3.4] with a closer attention to constants and the limit e —> 0 already taken. 

Lemma 3.10. Let x £ R n , r, <5 > 0, 0 <t<1 — \/y/2. Then 

limsup/r e ({u E > 1 — r} D B(x,r )) < 4 r p(B(x,r + <5)). 

£—>-0 

For all x £ R n , there are only countably many radii r > 0 such that p(dB(x,r )) > 0. This 
follows from the fact that p, is finite and that there are at most finitely many radii such that 
p(dB(x,r)) > 1/k, so that the union of those sets is countable. Thus for any r > 0 there is 
t £ (0, r) such that p(dB(x,t)) = 0. Letting S —> 0 at such a radius t (and using that the 
discrepancy measures go to zero) gives us the following result (compare also [DMR141 Lemma 9]). 

Corollary 3.11. For all x £ spt(p), r > 0 and t < 1/8 we have 

lirninf - C n ({|it e | < 1 — r} (~1 B(x , r)}) > 0. 

The following arguments rely more on the rectifiable structure of the measure p that we are 
approximating. Specifically, we introduce the diffuse normal direction by 

Vu e 

Ve:= W^\ 

when Vu £ 0 and 0 else. To work with varifolds, we introduce the Grassmannian G(n,n — 1) 
of n — 1-dimensional subspaces of R". We refer readers unfamiliar with varifolds or countably 
rectifiable sets to the excellent source |Sim83| : an introduction with different focus which is easier 
to find and covers most results relevant for us is [KP081 . Recall the following result. 

Lemma 3.12. |RS06l Propositions 4.1, 5.1] Define the n — l-varifold V e := p e (S> v e by 
V e (f)= f f(x,(^) ± )dp e V / £ C c (® n x G(n, n — 1)). 

JR n x G(n,n—1) 
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Then there is an integral varifold V such that V e —>• V weakly as Radon measures on R n x G(n, n—1) 
(varifold convergence). The limit satisfies 

Hv = fx, Hfi M < a 

where fxy the mass measure of V and H M denotes the generalised mean curvature of fx. In 
particular, W(fx) < d. 

The following result is a suitably adapted version of [RS06i Proposition 5.5] for our purposes. 
It shows that given small discrepancy measures and small oscillation of the gradient, a bounded 
energy sequence looks very much like an optimal interface in small balls. Using our improved 
bounds from Lemma |3.1| we can drop most of their technical assumptions. 

Lemma 3.13. Let 5, t > 0 and denote v e , n = ( v e , e n ). Then there exist 0 < L < oo depending on 
5 and r only and 7 > 0 depending on a, 5 and r such that the following holds for all x £ R". If 

(1) |u e (x)| < 1 — r and 

(2) ke\{B{x,4Le)) + f B(xALe) 1 - <7(4 Ley- 1 

then also the following two properties hold: 

• The blow up u e (y ) = u E (x + ey) is C 0ll / 4 -close to an optimal profile q on B( 0,3 L): 

| ± U e - q(y n — t \) |o,l/4,H(0,3L) < 8. 

The optimal profile q is the function from the lim. swp-construction and t\ = q~ 1 (u e (x)). 

• \u E (x,x n + i)| > 1 — r/2 for all Le < |t| < 2>Le, where x = (xi,..., x n _i) and u changes 
sign in between. 

Proof. Without loss of generality, we may assume that x = 0 and write B r := B{ 0, r). Recall that 
q'(t) = \/2W(q(t)) and lim^-too q(t) = ±1. Thus we can pick L > 0 such that \q(t)\ > 1 — r/4 
for all t > L. 

Assume for a contradiction that there is no constant 7 > 0 such that the results of the Lemma 
hold. Then for y - 7 —► 0, there must be a sequence u{ such that |u^(0)| < 1 — r, yV e {u{) < d + 1 and 


I ZiKBiL*) + [ 1 - «J 2 dtp < T j (4 Le) n -\ 

J B^Le 


but the conclusions of the Lemma do not hold. Considering the blow ups ft- 7 : B± b —* R with 
u 3 (y) = u J e (Ey) we obtain 


like in Lemma 3.1 


Since W 2 ’ 2 embeds 



||£t J '||2,2,B3i < Ca,n,L 

Hence there is u £ W 2,2 (B 3L ) such that 

w 7 —^ u in W 2,2 (B 3 l)- 
compactly into W 1,2 and L 4 , we see that 

I |Vm| 2 /2 — W(u) I dx = lim / I |Vu- J | 2 /2 — W(v?) I dx 

Jb 3L 

< lim £ 1 ^”|^|(R 4 l e ) 

j->oo 

< lim inf (4L)"-V 

j~* OO 


= 0 
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and when we set Vu = ( d\u ,..., d n -\u), we get 


>B S 


|V'S| dx = lim / |ViF|da; 
i^°° Jb 3L 


lim f \/\S7uj\ 2 — \d n ui\ 2 dx 
! ~>°° Jb 3L v 


< lim inf 

/ | VfPL 

Jl-fe) 

j-y °° . 

' BiL 

V V ) 


< 


lim inf (u n (4 L) n ) 1/2 ( f |VfP | 2 (l - {v° n ) 
^°° \Jb 4L 

1 ^ K) 2 d A 


< lim inf \J%Lui n ( e 1 n 

j-> oo \ 

< lim inf \J 8 L u> n yf 

j-yo° 

= 0 . 


' B 4Le 



1 

2 


Thus we can see that 

|Vft | 2 =2W{u) ) Vtt = (0,...,0 ,d n u). 

Clearly, this means that u(y) = p(y n ) for a function p with p' = ± \J2 W(p). Using that |m(0 )| < 
1 — t and the Picard-Lindeloff theorem on the uniqueness of the solutions to ODEs, we see that 
p(y n ) = ±<?(yn — y) for some y £ M. which can easily be fixed by the initial condition for p( 0 ). 

Since weak TU 2,2 -convergence implies strong C 0,1 ^-convergence in n = 2,3 dimensions, we see 
that there is j £ N such that the claim of the Lemma holds for u 3 e contradicting our assumption. 
Thus the Lemma is proven. □ 


To deal with the rectifiable sets in the next section more easily we prove a structure result for 
rectifiable sets. The result seems standard, but we have been unable to find a reference for it. As 
usual, we call a function on a closed set differentiable if it admits a differentiable extension to a 
larger open set. 


Lemma 3.14. Let M be a countably k-rectifiable set in R n . Denote by B the closed unit ball in k 
dimensions. Then there exist injective C 1 -functions R n with V/i ^ 0 on B such that 



and such that fi(B) fi fj(B) = 0 for all i j. 


Proof. According to [KP081 Lemma 5.4.2] or |Sim83l Lemma 11.1] there is a countable collection 
of C 1 -maps gi : R fc —> R" such that 


M C N U |J 9i (R fc ) 
1=1 


where TL k (N) = 0. Without loss of generality, N is assumed to be disjoint from the other sets. 
First we need to make the individual maps g, one-to-one. To do that, we define the set where 
injectivity fails in a bad way: 

Ai := {cc € R fe | V r > 0 3 y S B(x, r) such that g t (x) = gi{y) } . 
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Due to the failure of local injectivity, we see that the Jacobian J gi [x) vanishes on Aj. Since gi is 
a C' 1 -function, the set Di := J” 1 ^) is closed and by the Morse-Sard Lemma iFedfifll 3.4.3] 

n k ( 9i (Di)) = o. 

Set Ui := R k \Di. Now as in }EG92i Chapter 1.5, Corollary 2] we can use Vitali’s covering theorem 
[EG921 Chapter 1.5, Theorem 1] to obtain a countable selection of closed balls B? such that /, is 
injective with non-vanishing gradient on Bf for all j £ N and 


C k \ (J Bjj = 0. 

Since the boundary of a fc-ball has Hausdorff dimension k — 1, we could equally well take open 
balls. Since C 1 -functions map sets of £ fc -measure zero to sets of % fc -measure zero, we have shown 
that we can write 

OO 

M c IV U (J <ji(B°) 

J=i 

where T-L k {N) = 0, & : B —> R" is one-to-one, C 1 , and has a non-vanishing gradient everywhere 
on the closed ball B. The functions g m are obtained by rescaling suitable restrictions of gi from 
Bf to the unit ball. Finally, we have to cut out the sets that get hit by more than one function 
g m . Inductively, we define 


U m :=B°\~g- 1 U 9i(B)j ■ 

Finally, we use Vitali’s Lemma again to pick collections of closed balls B l m such that 

C k =0. 


1=1 


Rescaling the restricted functions from these balls and translating to the unit ball gives us the 
result. □ 


The proof of the following Lemma resembles that of the integrality of g in [RS061 Lemma 4.2]. 
It faces different challenges: while we do not need to prove multi-layeredness, we cannot zoom in 
on the tangent space since we need a macroscopic measure contribution to the double integral. 
Thus we need Lemma |3.14| to approximate macroscopically the structure of g. 

Lemma 3.15. Let (j) £ C°(1R.) such that <f> > 0 and <j)(u) dn > 0. If x £ spt (g), then 

liminf - f cf>(u e ) dx > 0 

e ^° £ J B (x,r) 

for all r > 0. 


Proof. Step 1. As usual, we assume that x = 0, g(dB(x,r/2)) = 0 and denote B = B(x,r/2). 
This means that all the e-balls of positive integral we are going to find will actually lie in B(x,r) 
and is a purely technical condition. Let £ be a small constant to be specified later. For further 
use, denote by B the closed unit ball in K n_1 . 

As g is an integral varifold, we know that spt(/r) is rectifiable. This means by Lemma 3.14 that 
there are countably many C 1 -functions fi : B — ► R" such that 


spt(M) C M 0 U (J ,h(B), H n ~\M 0 ) = 0 MB) n fj(B) = 0 
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for i ^ j- Since /z has second integrable mean curvature H 2 • /z < a, we can further use the Li-Yau 
inequality (3.3) to bound the maximum multiplicity of /z uniformly by 

a(n) 


< 


16 7r 


at least T~L n 1 -almost everywhere. Now since H n 1 (spt(/r)) < +oo we can find N gN such that 


7T- 1 f (spt (n) n B) \ (J MB)] < 

V i~ 1 / '-'max 


Since fi is injective and has non-vanishing tangent maps everywhere, M := Uf=i fi is a 

C'-manifold. We observe that 

K n ~\ spt(/i) n B\M)< 

”max 

and hence 

H(B\M)< C. 

Since the maps in question are smooth and the unit discs are orientable, for every i we can pick 
a continuous unit normal field to fi(B) (e.g. using cross products). Since the discs are compact 
and disjoint (thus a positive distance apart), the fields defined on each disc separately induce a 
continuous unit vector field on the union of their closures. 

Now we use the Tietze-Urysohn extension theorem to obtain a vector field X on B such that 
X = I'm on M and projecting on the unit ball we ensure |X| < 1. After an easy modification, we 
may assume that |X| = 1 on a neighbourhood of M. We then define 

G : R n x G(n, n — 1) —> R, G(x, S) = {X x ,u s } 2 


where vs is one of the unit normals to S. Note that G is continuous since X is. Using the non¬ 
negativity of G and the fact that Tj,/r = T X M for H n_1 -almost every x £ M D spt (fi) we interpret 
fi as dual to C 0 (R n x G(n,n — 1)) and observe 


(fi,G)= f 9{x)G(x,T x n)XH n 


L 

L 


6(x) G(x, T X M) d'H” -1 

f spt(/x)nM 

0(a:)d If- 1 

f spt(/x)nM 

= KM) 

> KB) - C- 

Step 2. By varifold convergence, we know that lim £ _ > o(/r e , G) = (/z, G) > n(B)—(, and |A|, \v e \ < 
1 so 

lirnsup / |l — (v e , X) 2 \ d/.i e = limsup / 1 — (v e , X) 2 dfi £ 

< limsup (h e {B) - (ne,G)) 

£—>-0 

<C- 

For 7, e, L > 0 we define the set 
U £ , 7 ,l := < x G B 


(4 Be) 


n— 1 


' B(x,4L/£) 


|l - {v e ,X) 2 I d n e > 7/4 > . 
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Let x±,... ,Xk be points in [/ Ej7i l being maximal for the property that the balls B(xi,4Le) are 
disjoint. Then by definition 


£ > [ 1 1 - {u e , X) 2 1 d Me > V [ \l-{v e , x ) 2 1 d/r E > K (ALe) n ~ l 7 / 4 . 

JB JB(xi,4Le) 

At the same time, we know that the balls B(xj, 8 Le) cover [/ Ej7i l because otherwise we could bring 
in more disjoint balls, therefore 


C n (U B ,y, L ) <K u n (8Le) n 

£ ~ £ 

< 4 £ u n (8Le) n 

~ 7(4L£) ra_1 £ 

= 2 n+i o Jn L(/'y. 


For a given 7, we choose £ = £(7) such that this is < n{B)/ 4. 

Step 3. Knowing that |£ e |(.B) —> 0, we can use the same argument as in the second step to 
show for 


V, 


e,7,L 


:= < x G B 


|6| (B{x,ALs)) 

(dLe)"- 1 


>7/2 


the estimate 

< MB)/4 

£ 


for all sufficiently small £ > 0. 

Step 4. Now choose U as a neighbourhood of M on which |Aj = 1 and r > 0 like in Corollary 
|3.11| satisfying 

liminf n s (U n {|u e | < 1 - r}) > 


This is easily achieved when i-i(M) > 3fj,(B)/4. Furthermore we take <5 <C 
small deviations of the optimal interface to behave similarly enough, L and 
and £ = £(7). Using steps one through three, we see that 


1 suitably small for 


7 as in Lemma 3.13 


liminf (jkl < 1 - t} n u \ (u E , 7 , L u u e , 7 ,l)) 

£—>-0 S 

> liminf £" (U e „, L ) _ C n (V e „,l) 
e ->0 e £ £ 

>3 /i(B)/4-/i(B)/4-MS)/4 

= ^(B)/4. 


Using the reverse argument of step 2, we can see that there are at least K points x\,... ,x k in 
{|%| <l-r}nU \ {U e ,-y,L U 14 i 7 i l) such that the balls B(xi,AL£) are disjoint with 


K > 


H(B) 

8 11+1 L n £ n ~ 1 ' 
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Step 5. To apply Lemma 3.13 we must “freeze” the coefficients of the vector field X to a single 
unit vector. We compute 


1 


(4 Le) 


n— 1 


' B(xi ,4:Le) 

l 


(l ~ ( u eiX) 2 ) — (l — (v e ,Xi) 2 ) d^ 


{ALe) n ~ 1 

1 


' B(xi , 4 Le) 


(u e ,Xi) 2 - {v E ,X) 2 dp, E 


(4 Le) 


n— 1 


{v e ,Xi-X) (u e ,X i + X)d^ e 

1 


f B(xi ,4Le) 

< |Xi + X\ c o(B{x it ALe)) ' |Xi - X\ C 0( B (xi,4Le)) U Le y-l ' 

< 2 C a ,L,n |Xi — X\co(B( Xi ,ALe)) 


’B(xi,4Ls) 


d[l E 


for all Xi such that \X{\ < 1. When we set Xi = X(xi), the last term converges to zero - so 
eventually it is smaller than 7/4 and 

1 f 


(4 Le) 


n— 1 


’ B{xi,ALe) 


1 - {v e ,XiY dn e < 7/2. 


Since Xi £ U, we finally see that |A'j| = 1 and Lemma 3.13 can be applied. 

Step 6. Since u e is C 0,1 / 4 -close to a one-dimensional optimal profile on B(xi,3Le) which 
transitions from —1 to 1, we see that for each s £ (—(1 — t),(1 — t)) there must be a point 
yi £ B{xi 1 3Le) such that u E (yi) = s. By Holder continuity, we deduce that 

[ dx > 9 E n 

J B(xi,3Le) 

for a constant 9 depending on the support of <j> and on d, n for the Holder constant. Since the balls 
are disjoint by construction, we can add this up to 

m „ 


1 r ^ \ 

- / <K“=) dx >~Yl 


J B(xi,3Le ) 


<j>(u e ) dx 


> -M9e r 
£ 


m( B)0 


~ 8" +1 L n 

> 0. 


This concludes the proof. 


□ 


We need one final result from geometric measure theory before we move on to our main results. 


Lemma 3.16. Let (X,d) be a metric space and K C X is compact. If K is not connected, then 
there exist two open sets U\,U 2 C X such that 

K C U\ U U 2 , K fl Ui Y 0 f or 2 = 1,2 and dist(f7i, U 2 ) > 0. 


Proof. Assume that spt(/u) is not connected. Then there exist relatively open non-empty sets 
Wi, W -2 C K such that 


K = W x U W 2 . 
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By definition of the subspace topology, W\ and W-j are also relatively closed. Since K is compact, 
they are even compact, so 5 := dist(M / i, W-f) > 0. We set U\ = {dist(-,Wi) < 5/3} and U 2 = 
{dist(-, W 2 ) < 5/3}. □ 

If K is the support of a measure /i, we can pick x t £ Wi for i = 1,2 since neither set is empty. 
As Ui is a neighbourhood of Xi, we see directly from the definition of the support of a measure 
that n(Ui) > 0 for i = 1, 2. 


3.4. Proof of the Main Results. Having dealt with the necessary auxiliary results, we can 
proceed to prove our main results. We use the terminology of Lemma |3.9[ 


Proof of Theorem 2.1 Compare also [NT07 1 Lemmas 3.1, 3.2 and Proposition 4.2], Let Q! <e 
K™ \ spt(//) and t > 0. Assume there is a sequence x E in fl' such that |u e (a; e )| < 1 - r or 
|u e (x e )| > 1 + r. By definition, x e £ Q and using compactness, there is x £ tt' n such that 
x e x. 

Let r > 0 such that B(x,3r) Cl" \ spt(//). Due to convergence, B{x s ,r) C B(x,2r) for all 
sufhciently small e > 0. If we use the rescaling property from Lemma |3.4| and the minimisation 
property of Lemma |3.9| for n = 2, we see that 

1 


a(B(x,2r))> - liminf (r a e {B(x e ,r)) + r 1 /j. £ (B(x s , r))) 
r e->0 

> - liminf inf TAi t) 

r e->o ueY 2 

> 0 o /r 


with i = ejr as in Lemma 3.4 Letting r —> 0, we obtain a contradiction. 


□ 


If n = 3, neither the rescaling property of Lemma 3.4 nor the minimisation property of Lemma 
3.9| hold in as strong formulations. As we sketched in Remark |2.8[ uniform convergence is, in fact, 
false. 


Proof of Theorem \2.S\ : We assume that u E —> u strongly in L 1 ^) and that /i e /r for our sequence 
or a suitable subsequence (not relabelled). Let D (-1, 1), then u e 1 (D) <e Q. We can consider 
u 7 1 (D) instead without changing the Hausdorff limit to conform with standard approaches. By 
the usual compactness results (see e.g. (KP081 Theorem 1.6.6]), there is a compact set K such that 
a further subsequence of uf 1 (D) converges to K in Hausdorff distance. K can be computed as the 
Kuratowski lower limit 

K = {x £ R" | 3 x E e m/ 1 (D) such that x E — ► x). 

Step 1. We first show that K C spt(/u) if n = 2. Assume x G K \ spt (fi). Then there exists 
r > 0 such that B(x,r ) <e R 2 \ spt(/i). This means already that |u e | —> 1 uniformly on B(x,r) 
showing that uf 1 {D) n B(x,r ) = 0 for small enough e, contradicting our assumption. 


Step 2. Now we prove that spt (/r) C K in n = 2,3 dimensions. It suffices to show that for 
x £ spt (/j,), s £ D and r > 0 we have w/ 1 (s) D B{x, r) ^ 0 for all sufficiently small e, which implies 
Hausdorff convergence. 



If n = 2, the uniqueness of the limit also shows convergence for the whole sequence. □ 


We will now proof our main statement about connectedness. 
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Proof of Theorem 2.5. The proof is structured as follows. First, we show that we can find neigh¬ 
bourhoods of connected components which have positive distance with respect to the usual metric 
on R". Then we need to show that they also have positive distance with respect to the pseudo¬ 
metric d F ( Ue f Intuitively, this makes sense since any connecting curve should have to leave the 
interfacial layer between the two sets. This is simple if n = 2 and slightly more technical if n = 3. 
We proceed in the spirit of Theorem 2.3 and show that u e converges to ±1 even on curves. 


Without loss of generality, we may assume that there are — 1 < 9\ < 0 2 < 1 such that {cf > 0} C 
(di, 9 2 ) and F > 1 outside (9 1 , 62 )- This is only a minor assumption and could easily be removed, 
but simplifies the proof. 

Step 1. Assume that spt (/i) is not connected. Since spt(n) is compact, according to Lemma 
|3.16| there are disjoint open sets U\, U 2 such that 


spt(/x) C U\ U U 2 , h{Ui) > 0, i = 1,2, 5 := dist(f7i, U 2 ) > 0. 


Now 


lirninf C F (u F ) > liminf 
e—>o v ' - e—>0 



<f>(uAx)) da; • liminf 

e->0 



4 >{u e (y))dy 


■ liminf dist F ^“^(J/i, U 2 ). 
£—>-0 


Since the first two factors are strictly positive according to Lemma 3.15 it suffices to show that 
liminf £ _> 0 dist F( “ e) ([/ 1 , U 2 ) > 0. 

Step 2. For a contradiction, assume that dist F ^“ e ^ (U\, U 2 ) —> 0. Pick a sequence c e such that 


c £ -> 0 but still 
such that 


dist F( " e) (!7i,H2) 


0. Then there exist a connected set K e and points x e ,y s € 
x E e K e n dU x , y e £K e ndU 2 , [ F(u e ) d'H 1 < c e . 

JK r 


If n = 2, we know that \u e \ —>• 1 uniformly on \ {U\ U U 2 ), so in particular u e ^ [9\,6 2 ] on 
K e C f2 \ (Ui U U 2 ) and 



F(u e ) d n 1 > U l (K e \ (t/i U U 2 )) > 6 > 0 


since K e connects U\ to U 2 . This is a contradiction to our assumption. In the case n = 3 we need 
a further argument. 

Step 3. In this step, we will use the competition between the distance function driving u e away 
from ±1 along K e and the energy bounds in three dimensions. 

Precisely, knowing that the F(u £ )-weighted length of K e is small, we will construct a set of 
N e 1 points Xi i£ on K e such that u e (x j )£ ) € [0i, 6 2 ] for all * = 1,... N and such that the balls of 
radius e 1 / 2 around these points are disjoint. The large number of balls implies that one of them 
has to satisfy the conditions of Lemma |3.9[ This allows us to show that a must have an atom of 
fixed minimal size close to this point. 

A more careful construction allows us to localise the argument on short segments of K e such 
that we may construct any finite number of atoms with the same fixed minimal size. This gives 
the desired contradiction. 

Take a subsequence realising the liminf. The 1-Lipschitz map 7r(x) := dist(x, U\) maps K e to 
a connected set containing 0 = 7t(x £ ) and 5 = 7 r(y e ), so [0, A] C 7r(A' £ ). Furthermore, 7r _1 (0, 5) C 
M 2 \ (U\ U U 2 ) since dist(C/i, U 2 ) = 5. Take the set 


K' e ■= {t £ [0, <5] : 3 x G K e such that t = ir(x) and u e (x) ^ [0i, 0 2 )} 
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of points whose pre-image contributes a lot to the weighted length of K e . Then 

H 1 (K' e )<H 1 (K e n{u e <t[e 1 ,0 2 ]}) 

< f F(u e ) AH 1 

JK r 


>K e 
< Cr. 


Pick M intervals 


h = 


2k — 1 k 
2 M ’ M d 


inside [0,<5]. Fix 1 < k < M. When e is so small that c E < we deduce that 


(3.5) 




Take a maximal collection ti j£ ,..., tiv e ,£ i n Ik with the property that the balls l?(fj jE , e 1 / 2 ) are 
mutually disjoint. Then N e > 16A / i / 2 — 1 since otherwise 


N e 


I k \ K’ e c U B{U,e, 2e 1/2 ) =► n\h \ K’ e ) < N e ■ Ae 1 ' 2 < 


i=1 


16-Me 1 / 2 


- 1 • 4^/2 


< 


4M' 


For each t ii£ we pick Xi >£ £ _?f £ fl 7 r _ 1 (t ij£ ) which satisfies by construction u e {x^ e ) £ [ 0 i,# 2 ]- 
Since the projection map is contracting, also the balls B(xi )E , e 1 / 2 ) are mutually disjoint. Fix any 
0 < 77 < 1/2. Now assume for a contradiction that all points Xj >£ satisfy 


a e (B(xi, e ,£ 1/2 )) + ( 
J E 


W\u E f 


Ax > e r ‘ 


for all* = 1,..., N E . Thus 


2a E (f l) > a E {fl) + f 
./( 


W’{u e ) 2 * 


-da; > N f £ v > 


{ I U £ I > 1} 


~ 3 £ 16M 


£ V -i/2 _ £V ~ £V -i/2 


due to (3.1 1 . This grows beyond any bound as e —> 0, so there must by 1 < i e < N e such that 

W'(u E ) 2 


Oi e (B(Xi t s, £ 1/2 )) + [ 
J E 


da; < £ v . 


s(x i , e , E 1 /2)n{|u e |>i} £ 


Using compactness and a subsequence, there exists x k £ 7 r - 1 (/fc) such that Xi ejE —>• x k . We can 
achieve this simultaneously for all k = 1,... ,M. The balls B(x k ,r ) are disjoint for 0 < r < s 
since the intervals Ik have distance S/2M and y{B(x k , r)) = 0 since 


4 M 


dist(xfc, spt(/z)) > dist(U \ TJ\ U U 2 , spt(//)) > S. 

Fix such r and consider e = e/(2r) and 

u £ : U(0,1)-)• R, u e (y) = u e Q • y + x i>£ J 

with u £ (0) £ [0i, # 2 ]- We use the terminology of Lemma |3.4[ The associated measures a E satisfy 

r W’{u E ) 2 


A e (B( 0,1)) 


Ax < £ n 


/s(o,i)n{|u e |>i} 
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so by Lemma pT~TT| 

a(B(x k ,r )) > liminf a E (B(x k , 3r/4)) 

£->0 

> liminf a e (B(xi E , r/2)) 

£—>■0 

= liminf [a E (B(x ie>E ,r/ 2)) + (r/2) 1- ”/x £ (-B(:Ei e , £ , r/2))] 

= lim inf (/i e + a E )(B(0, 1)) 

£—>•0 

> 00- 


Taking r —> 0, we see that a({x k }) > 9 0 . So in total, a(Q) > M9 q for all M £ N where 9 0 depends 
only on 9\ and 02- Thus we have reached the desired contradiction. □ 


Now Theorem 2.6 is an obvious consequence of Theorem |2. 5 


Proof of Theorem \2.9\ Let u e be a sequence such that £ e (u e ) is bounded. Then in particular 
|^ e (n) — S\ < e' 7 / 2 , so fi E (M. n ) = /n e (fi) is bounded and /r e —*■ /r for some Radon measure /Lt for 
this and other properties see [EG921 Chapter 1], Clearly 


/,i(f2) > limsup/Lt e (n) = S 
£—^0 


and on the other hand 

li(n) < n(R n ) < lim inf /i e (R n ) = S 
£->-0 


so /x(f2) = S. If U = M ra \ 12, we have 


/x(17) < liminf /i e (17) = 0, 
£—>■0 


so spt(/i) = Her open „(u)=o u ° Since £ E (u E ) is bounded, we have C e (u E ) —> 0, so due to 

Theorem |2.5| spt(//) is connected. □ 

Finally, the same arguments give convergence H 1 -a.e. on fixed curves. 


Proof of Theorem \2.3[ This proof is a simpler version of that of Theorem 
everywhere on K , there is a set 


2.5 


If \u e \ -ft 1 "H 1 -almost 


K' s := {x € K | lim sup | |u e | — 1 I > s} 
£—>-0 


with 7 ft{K' s ) > 0 for some small s > 0. We do not need the projection onto the first coordinate 
and can directly define x i e in the fixed set K e = K' s without passing through t ii£ , the rest works 
analogously. □ 


4. Computer Implementation 

In a computer simulation, we try to find local minimisers of the £-problem in n = 2 dimensions. 
In the implementation, we follow a finite element approximation of the time-normalised L 2 -gradient 
flow of £ e with e = 1.5- ICC 2 , n = 1, a = 2. The pictures below were obtained on a triangular mesh 
with elements of diameters approximately h = 6 • ICC 3 and approximately 250.000 i7 2 -conforming 
subdivision surface basis functions supported in the two-ring around an element. Time-stepping 
was semi-implicit with a time-step of r = £ • 10 -5 . The distance function d F ( Ue ' ) and its gradient 
are implemented via Dijkstra’s algorithm in a fashion similar to the one of |BCPS10] , The initial 
condition is the same for all simulations and can be seen in Figure [l] on the left; the domain is the 
disc of radius 1. There was no penalisation of the discrepancy measure in our simulations. 

For practical purposes, we use two functions with support close to 1 and — 1, respectively, 

rather than just one 4>. Quite obviously, our proofs easily extend to that situation. By keeping 
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level sets close to the edges connected, we create barriers that prevent the interface from splitting 
apart early in the process. The implementation will be described in greater detail in a forth-coming 
article l )\\ Hi . 

We see in Figure [l] that without the inclusion of the topological term, the transition layer 
disintegrates into several connected components along the gradient flow of W £ + e~ a ( S E — S) 2 . 



Figure 1. Gradient flow of W £ +e a (S e — S) 2 . From left to right: phase field u 
for approximately t = 7.5 • 10 —5 , t = 3 • 10 -4 , t = 7.5 • 10 -4 and t = 1.8 • 10 -3 . 

To compare implementations of topological side conditions, we include the topological term 
suggested in j DMRllj . which penalises a deviation of a diffuse signed curvature integral from 2i r 
in the simulation. This term prevents the initial pinch-off, but at a later time, the interface will 
pinch off in a more complicated way which keeps the diffuse winding number close to 27r. The 
phenomenon is a simultaneous pinch off at several points as seen in Figure [2] The far right plot in 
Figure [2] illustrates the diffuse curvature density as distributed along the curve at pinch off time. 
We can observe the formation of a circle with negative total curvature ss — 2ir (due to the phase 
field switching in the other direction from +1 to —1), and two components with total curvature 
ss 27 t so that the total curvature of the whole interface stays close to 2 ir. 



Figure 2. Gradient flow with penalty on a diffuse winding number as suggested 
in |DMRll| . From left to right: phase field u for approximately t = 3 • 10 -4 , 
t = 7.5 ■ 10 -4 and t = 1.8 • 10 -3 , then a plot of the diffuse winding number density 
denoted T at time t = 1.8 • 10 -3 . 

In Figure [3j a flow for 8 e with the additional term of C e on the other hand can be seen to stably 
flow past those singular situations. 

Comparing the three scenarios above, we observe that there is virtually no difference in the 
plots at time 3 ■ 10 -4 and that the plots for both modified (penalised using either the old or the 
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Figure 3. Evolution including our new topological penalty term C e . Top line, 
from left to right: phase field u for approximately t = 3 • 10 -4 , t = 7.5 • 10 -4 
and t = 1.8 • 10 -3 , then a plot of the diffuse Willmore energy density (denoted 
W here) of the initial condition, bottom line, left to right: phase field u and 
diffuse Willmore energy density first for approximately t = 6.6 • 10“ 3 and then for 
approximately t = 3.6 • ICC 2 . 


new method) functionals at time 7.5 • 10 -4 still look very similar. It can thus be argued that the 
topological condition does not affect the shape of the curve in a major way except when it has to 
in order to prevent loss of connectedness. 

In Figure [3j we see non-trivial geometric changes along the gradient flow for later times. This 
demonstrates the necessity of continuing the flow beyond the critical times. 

It should be emphasised that our focus is not on implementing a scheme to approximate Willmore 
flow using phase fields but on finding minimisers of the diffuse interface problem using a gradient 
flow. Existence of Willmore flow for long time and topological changes along it are still an open 
field of research. 


5. Conclusions 

In this paper, we have developed a strategy to enforce connectedness of diffuse interfaces. The 
strategy fares well in applications and can efficiently be implemented and seems to be more generally 
applicable to a wider class of problems. We claim that our results can be extended to the following 
situations. 

• We can include a soft volume constraint like 

ir G/,"* +idi ) 

for continuous functions F > 0. 
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• Another popular constraint compatible with our functional and results is minimising a 
distance from a given configuration as 


A e (u) = / \u — g\ dA 
Jn 

where A is a finite Radon measure on 11 and g £ L 1 (fl). This functional originates in 
problems in image segmentation, but in our context it can be understood as prescribing 
certain points to lie inside or outside the membrane according to experimental data. 

• Using [BM10I . Theorem 4.1], we could take Bellcttini and Mugnai’s approximation of the 
Helfrich energy 



for x G (—2,0) in place of the diffuse Willmore energy W £ . Here v u , e is the usual Willmore 
density associated with u and v u = Vu/|Vzt| is the diffuse normal like in Lemma 3.12 


We can use the same modelling techniques for a finite collection of membranes given by 
ul,.... u^ inside an elastic container given by U e . The governing energy could be composed 
of a sum of the individual elastic energies S E and interaction energies I E like 

1 


= 


/« 

Jn 


+ l) 2 (u» + l) 2 dx 


which prevent penetration of the phases u l e ss 1 and u 3 e ~ 1 or, in a slight variation, enforce 
confinement of u\ « 1 to U e ss 1. 
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